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Modification of Born impurity scattering near the surface of d-wave superconductors 

and influence of external magnetic field 
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We study the influence of self-consistent Born impurity scattering on the zero-energy Andreev 
bound states near the surface of a d-wave superconductor with and without an externally applied 
magnetic field. Without an external magnetic field we show that the effect of Born impurity scat- 
tering is stronger at the surface than in the bulk. In the presence of an external magnetic field the 
splitting of the zero-energy Andreev bound states is shown to have a nonmonotonous temperature 
dependence. Born impurity scattering does not wash out the peak splitting, but instead the peak 
splitting is shown to be quite robust against impurities. We also show that a nonzero renormalization 
of the pair potential appears near the surface. 
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I. INTRODUCTION 

At the surface of d-wave superconductors zero-energy 
Andreev bound states may appear depending on the 
orientation of the d-wave with respect to the surface 
normalii^i^i^. Experimentally, these states can be ob- 
served as zero-bias conductance peaks in the tunneling 
conductance^iii^. It is well known that surface rough- 
ness, surface disorder^, or diffuse scattering at the sur- 
face leads to a broadening of these states. Also, impurity 
scattering in the bulk of the superconductor is known 
to broaden the Andreev bound statesiSiii. In the pres- 
ence of an external magnetic field the screening current 
leads to a splitting of the Andreev bound state o^^i^^'^'^ . 
In this case, a counter-flowing paramagnetic current is 
generated by the surface states, which increases with de- 
creasing temperature and may even lead to a reversal of 
the current flow at the surface resulting in an anomalous 
Meissner effecti^ii^. This effect has recently been shown 
to have a strong influence on the Bean-Livingston surface 
barrier for entrance of vortices into the superconductor^^. 

In the present work we investigate the influence of bulk 
impurity scattering on the broadening of the surface An- 
dreev bound states and the splitting in an external mag- 
netic field. It has been shown previously that impurity 
scattering in the Born limit is much more effective in 
broadening the Andreev bound states than impurity scat- 
tering in the unitarity limi t^°'^^ . In the high-Tc cuprate 
compounds it is believed that scatterers within the Cu02 
planes act as unitary scatterers and thus should have lit- 
tle influence on the Andreev bound states. However, it 
has been recognized recently that scatterers sitting be- 
tween the Cu02 planes are poorly screened and act as 
Born scattereraiiiSiii^. These impurities are thus ex- 
pected to have a dominating influence on the broadening 
of the Andreev bound states. For these reasons in the 
present work we will focus on the influence of impurity 
scattering in the self-consistent Born approximation. We 
will show that in this case impurity scattering around 
zero energy is significantly increased near the surface as 



compared to the bulk, leading to a larger broadening of 
the Andreev bound states than expected from the scat- 
tering rate in the bulk. The situation changes completely 
in the presence of an external magnetic field, however. 
The splitting of the Andreev bound states turns out to 
be quite robust against Born impurity scattering. 

In the bulk of a d-wave superconductor the renormal- 
ization of the pair potential due to impurity scattering is 
known to disappear. However, this is not generally the 
case near a surface because of broken translational sym- 
metry. Here, we will show that a nonzero renormalization 
of the pair potential appears near the surface unless the 
orientation of the surface is highly symmetric with re- 
spect to the orientation of the d-wave. Also, in the pres- 
ence of an external magnetic field the renormalization of 
the pair potential becomes nonzero. 

In Section [11] wc will describe our numerical approach. 
In Section Hill we will first study impurity scattering near 
a surface of a d-wave superconductor in the absence of 
an external magnetic field. Section IIVI presents results 
for a superconductor without impurities in the presence 
of an external magnetic field, and Section |V] deals with 
both impurity scattering and the presence of an external 
magnetic field. 



II. 



NUMERICAL APPROACH 



The geometry under investigation is shown in Fig. [T] 
in the halfspace x > 0, we assume to have a supercon- 
ducting area of d-wave type. For a; < 0, an external 
magnetic field B = Be^ is applied parallel to the z-axis. 
For simplicity, we consider a cylindrical Fermi surface 
with the c-axis oriented parallel to the z-axis. Also, we 
assume that the external magnetic field remains smaller 
than the field of first vortex penetration. Therefore, we 
can assume translational invariance both along the y-axis 
as well as along the z-direction. The relation between the 
current density in the superconductor and the vector po- 
tential in Coulomb gauge is given by Maxwell's equation: 




FIG. 1; (Color online) Sketch of the investigated geometry, 
showing the pair potential and the direction of the vector 
potential relative to the surface of the superconductor. The 
angle a determines the relative orientation of the d-wave with 
respect to the surface normal. 



Here, the renormalized Matsubara energies and pair po- 
tential are given by: 

ie„(i?(s),e„,T) = i£„ + -z;j.-A(i?(s))-S«(i?(s),e„,T) 

c 

A(i?(s),e„,T) = A(i?(s),r)cos(2(0-a)) + 
E^(^(s),e„,T) 

where E*^ and E^ denote the diagonal and ofF-diagonal 
self energies due to impurity scattering. The pair poten- 
tial is given by: 

A(r, T) = VN^^T ^ (cos(2(0 - a))/(f, kp, ze„))Fs(5) 

The brackets {■ ■ ■)fs denote an angular average over the 
cylindrical Fermi-surface. By solving the Riccati equa- 
tions along real space trajectories R{s) running paral- 
lel to the Fermi velocity vp^ the normal and anomalous 
propagators are found from: 



AA = —3 
c 



(1) 



In this gauge the vector potential is directly proportional 
to the superfluid velocity Vs- The boundary conditions 
for this second order differential equation are determined 
by the behavior of the magnetic field: it penetrates into 
the superconductor continuously and decays to zero in 
the bulk. 

Our calculations are based on the Eilenberger 
equations22i2i. These equations can be transformed into 
Riccati type differential equations for two scalar com- 
plex quantities a{s) and h{s) along real space trajectories 



22. 



SVp 



hvF^a{s) + [2l,^{s) + K\s)a{s)\a{s)-K{s) = 
hvF^Ks) - [2?„(s) + A{s)b{s)]b{s) + At(s) = 0(2) 

OS 

Here, vp is the Fermi velocity and vp the unit vector 
in the direction of the Fermi velocity. The initial values 
for solving the Riccati equations for a{s) and b{s) are 
obtained from the fact, that their variation vanishes in 
the bulk: 



a(— oo) 



A(-oo) 



6l + |A(-oo)P 



(3) 



H^»)-H)-^^44g4> fim)-T,^^i^) 



1 + a{s)b{s) 



1 + a{s)b{s) 



From the propagators, we can instantly derive the cur- 
rent density j(r) and the local density of states (LDOS) 

N{E,f): 



N{E,r) 
No 



-lm{g{r,kFjen^E + iO+))FS (7) 



where Nq is the normal state density of states at the 
Fermi level, and 



j{r,T) = 2eNoVFkB'n'T ^ {vf ■ g{r,kF,ien))FS (8) 

\en\<^c 

The zero temperature London penetration depth Al is 
given by the expression A^^ = ^e^NoVp. Throughout 
the work we will quote A^ relative to the zero tempera- 
ture coherence length without impurities ^o = 7a7§^oT' 
i.e. the parameter k = Xl/£,o- Magnetic fields will be 
given in units of the zero temperature upper critical field 
Bc2, which for a bulk d-wave superconductor is given by 
Bc2 = 0.49 15^^ for a cylindrical Fermi surface. 



In our model, we include the effect of impurity scat- 
tering in the self-consistent Born approximation. As 
pointed out above, Born impurity scattering is expected 
to cause a stronger effect than scattering in the unitar- 
ity limii^iSiii. For simplicity, we will restrict ourselves to 
s-wave scattering. In this case the impurity self energies 
are given by: 



b{+oo) = 



At(+oo) 



|A(+oo)|2 



(4) 



E^(f,e„,r) = _(/(r,fcj.,e„))f5 (9) 

S<=(r,e„,T) = l-{g{f,kF,e„))FS (10) 




FIG. 2: (Color online) Bulk value of the pair potential as a 
function of mean free path l/^o for different temperatures. 




FIG. 3: (Color online) Negative imaginary part of self energy 
E*^ as a function of energy E for different distances from the 
surface at an orientation angle of a = ^ . The temperature is 
T = O.lTc and the mean free path I — 2.7^0- 



where t is the scattering lifetime in the bulk and is given 
by: 



III. IMPURITY SCATTERING NEAR A 
SURFACE OF A d-WAVE SUPERCONDUCTOR 



':7TNon,\Vo\ 



where rii is the impurity concentration, Vq is the strength 
of the impurity potential. Throughout this work, we will 
quantify the impurity scattering in terms of its mean free 
path I — vpT relative to the zero temperature bulk co- 
herence length in the superclean limit ^q- Iii the bulk 
of a d-wave superconductor the angular average of the 
anomalous Green's function / over the Fermi surface van- 
ishes, because positive and negative contributions cancel 
exactly. This leads to a vanishing renormalization E^ 
of the pair potential. However, as we will show below, 
this does not generally hold anymore in the vicinity of a 
surface. 

Numerically, we start with estimated functions A(a;) 
and A{x). These are used to solve the Riccati equations 
^ along all real space trajectories with specular reflec- 
tion on the surface a; = 0. From the solutions we find the 
propagators Eq. ([5]). These are used to obtain the self 
energies Eq. ^ and (jlOp. the current density Eq. ([5]), 
and the updated pair potential Eq. ([5]). Finally, integra- 
tion of Eq. Il]) yields an updated vector potential. This 
procedure is iterated until the functions A(a;) and A{x) 
converge. Note, that the self energies E*^ and E^ and 
the propagators g and / are calculated self-consistently 
this way. After convergence, a final iteration is run, in 
which all equations are solved directly for real frequen- 
cies ?e„ -^ E -\- iQ^ in order to perform an analytic con- 
tinuation for the local density of states and the self en- 
ergies. Only for the calculations without impurities we 
have added a small imaginary part of Q.^QlksTc in order 
to regularize the poles of the propagators. 



In this section we consider a superconductor with dif- 
ferent impurity concentrations for the case that no ex- 
ternal magnetic field is applied. It is well known that 
impurities in bulk d-wave superconductors lead to pair 
breaking, which implies a decrease of the bulk order pa- 
rameter. Fig.[2]shows the bulk value of the pair potential 
as a function of the mean free path I for different temper- 
atures. When the mean free path becomes comparable to 
the finite temperature coherence length, the pair poten- 
tial vanishes. Consequently, for d-wave superconductors, 
there exists no real dirty limit. 

Near the surface Andreev bound states are absent for 
d-wave orientation a = 0. When the angle a is increased, 
the spectral weight of the Andreev bound states gradu- 
ally increases until it reaches a maximum at a = 7r/4. In 
the following we will first focus on an intermediate angle 
of a = 7r/8, where the spectral weight is neither absent 
nor fully developed. 

The local quasiparticle scattering rate is given by the 
negative imaginary part of the normal self energy -Im S . 
In Fig. [3] we show the energy dependence of -Im Tp 
for different distances from the surface at a temperature 
T ~ O.lTc and a mean free path of I — 2.7^o- From this 
it can be seen that there is a significant variation of the 
quasiparticle scattering rate at the Fermi level _E = as 
a function of the distance from the surface. For this set 
of parameters, at the surface the quasiparticle scattering 
rate is about 12 times larger than in the bulk. Physically, 
this effect can be understood from Eq. PH]) . At the sur- 
face in the presence of the Andreev bound states there 
is a larger phase space of low energy states available for 
scattering. This makes impurity scattering more effective 
at the surface than in the bulk. 

In Fig. m the impurity dependence of the local den- 
sity of states at the surface is shown for the same an- 
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FIG. 4: (Color online) Local density of states at the surface 
for different impurity concentrations. The temperature is T = 
O.lTc and the orientation angle a = ^. 
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FIG. 6: Zero energy density of states at the surface as a func- 
tion of \j J- for two different orientations a = ^ (solid) and 
a = f (dashed). Temperature is T = O.lTc. 
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FIG. 5: The local density of states at the surface for temper- 
ature T — O.lTc, orientation angle ce ~ -^ and mean free path 
I = 2.7^0. The dashed curve shows the LDOS when the bulk 
values of the self energy are used. The solid curve, on the 
other hand, shows the LDOS when using the self consistent 
solution for the self energy. 



gle a and temperature. Increasing the impurity con- 
centration results in a decrease of the height of the 
zero energy peak and a broadening of its width. The 
peaks seen near lAksT^ and O.O/csTc in the absence of 
impurity scattering can be interpreted as follows: the 
peaks near ±.\AkBTc are related to the bulk gap times 
cos 2a. They are coming from quasiparticles, which ap- 
proach the surface perpendicular, as a look at the mo- 
mentum resolved data shows. In contrast, the peaks near 
zbO.QfcsTc are caused by grazing angle quasiparticles. In 
our self-consistent calculation the local gap near the sur- 
face is much smaller than in the bulk. The grazing angle 
quasiparticles mostly experience the reduced surface gap 
value, creating a gap edge around COfc^T'c- In the pres- 
ence of impurities these gap features are quickly washed 
out. 



In order to illustrate the change of the local density of 
states at the surface due to impurity scattering, in Fig. [5] 
we compare the local density of states at the surface for 
/ — 2.7^0 with a hypothetical calculation, in which we 
have used the bulk value of Yp at the surface. Clearly, 
the zero energy peak is much sharper when the bulk Yp 
is used for the calculation. 

These results show that the influence of Born impurity 
scattering is much stronger at the surface than in the bulk 
due to the presence of the Andreev bound states. Their 
presence creates a larger number of available scattering 
channels, which in turn leads to a stronger broadening 
of the Andreev bound states. This self-consistent broad- 
ening can be illustrated by looking at the peak height of 
the local density of states at zero energy N{E — Q)/Nq. 
On the one hand the peak height scales approximately 
with the inverse of the local quasiparticle scattering rate: 



N{E = 0) 



kuTr 



-ImI]G(£; = 0) 



On the other hand the local quasiparticle scattering rate 
is determined by the peak height via Eq. PU)) : 



lmS^{E = 0) = 



1 N{E^O) 



2r TVo 
Solving for the peak height leads to the expression 

N{E = 0) 



This result is in good agreement with the numerical result 
shown in Fig.[6l It shows that the scaling behavior of the 
peak height is ^ vl instead of the ^ I behavior one would 
have expected from bulk scattering, leading to a stronger 
impurity effect near the surface. 




FIG. 7; (Color online) Real part of self energy E^ at the 
surface vs. orientation angle a for different energies. The 
temperature is T = O.STc and the mean free path I — 2.7^o- 
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FIG. 8: (Color online) Real part of the self energy E^ vs. 
distance for two different energies E. The inset shows the 
energy dependence of this self energy for different distances 
from the surface: a;/^o = (solid line), a;/5o = 5 (dashed 
line), and x/^o = 20 (dash-dotted line). In both cases, the 
orientation is given by a = ^, the mean free path I = 2.7^o, 
and the temperature is T = O.lTc. 
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FIG. 9: (Color online) Current density distribution for ori- 
entation a = J at temperatures T = 0.1 Tc (solid), 0.4 Tc 
(dashed), 0.7Tc (dash-dotted), 0.9 Tc (dotted). The exter- 
nal magnetic field is B^xt = 0.02i3c2. The inset shows the 
temperature dependence of the surface current density for 
orientation a — j and the same value of the external mag- 
netic field. Here, the current density has been normalized to 

Jext — 4^ ^ ■ 



of S^ for d-wave superconductors is not generally true 
anymore near the surface. Near the surface translational 
invariance is broken, which makes trajectories with dif- 
ferent momenta kp inequivalent, because they experience 
different pair potential landscapes. Except for special ori- 
entations a of the d-wave with respect to the surface this 
leads to finite values of the anomalous self energy E^. In 
Fig. [7] we show the real part of S^ at the surface as a 
function of the orientation angle a for different energies. 
It can be seen that E^ vanishes for integer multiples of 
7r/4. Fig. [5] shows how I]^ varies with the distance from 
the surface and energy (inset) for a = tt/8 decreasing to 
zero in the bulk. 



IV. INFLUENCE OF AN EXTERNAL 
MAGNETIC FIELD IN THE CLEAN LIMIT 



In an isotropic s-wave superconductor Anderson's the- 
orem asserts that the renormalization of the pair poten- 
tial due to the anomalous self energy E^ exactly com- 
pensates the renormalization due to the normal self en- 
ergy T,'^ , such that the density of states and Tc remain 
unaffected by impurity scattering. This does not nec- 
essarily hold anymore in an anisotropic superconductor, 
however—'^. In the bulk of a d-wave superconductor the 
anomalous self energy E^ is known to vanish. This is 
clear from Eq. ©, because the Fermi surface average 
leads to cancellation due to the sign change of the d- 
wave. Ultimately, this is the reason why nonmagnetic 
impurity scattering is much more destructive to uncon- 
ventional superconductors than to conventional ones. It 
has not been noted before, however, that this vanishing 



In this section we will discuss the influence of an exter- 
nal magnetic field in the absence of impurity scattering. 
In particular we focus on the case a = j^ where the 
spectral weight of the Andreev bound states is strongest 
and K = 10. This value of k is modest in comparison 
with K values of hole doped high- Tc cuprates, but may 
be relevant for some low Tc electron doped cuprates2i. 
The influence of the anomalous Meissner effect becomes 
more pronounced for small values of k and here we wish 
to illustrate a peculiar effect that occurs in this range of 
parameters. 

In the following we have set the external magnetic field 
to Bcxt = 0.02Bc2- Fig. [HI shows the current density 
for selected temperatures as a function of the distance 
from the surface. In a distance up to 3^o from the sur- 
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FIG. 10: (Color online) Magnetic field as a lunction of the 
distance from the surface for orientation a = j and different 
temperatures. The external magnetic field is Bext = 0.02_Bc2. 



face (which is of the order of the spatial extension of 
the bound states), the current is flowing opposite to the 
screening currentA^'ii. This anomalous Meissner current 
persists throughout the full temperature range between 
O.OlTc and 0.9Tc, as demonstrated in the inset of Fig.[9l 
While it is nearly vanishing for temperatures close to Tc, 
it saturates near zero temperature . 

The magnetic field distributions resulting from the cur- 
rent distributions are shown in Fig. 1101 With the anoma- 
lous Meissner current flowing, the magnetic fleld initially 
increases before the normal Meissner screening sets in 
and eventually screens out the magnetic fleld exponen- 
tially. This initial increase occurs again up to a distance 
of ^ 3^0 from the surface. With the value of k = 10 we 
have used here, the fleld increases by more than a factor 
of two relative to the external fleld. Qualitatively it is 
clear that this fleld increase becomes more pronounced 
for smaller values of k, because a smaller penetration 
depth results in larger current densities, as seen from 
Eq. (0). 

Fig. [TT] shows the modulus of the vector potential at the 
surface as a function of temperature. It can be seen that 
the temperature dependence is nonmonotonous. The vec- 
tor potential increases both towards low temperatures as 
well as towards Tc- The behavior near Tc is easily un- 
derstood from the temperature dependence of the pene- 
tration depth, which diverges near Tc- Since the vector 
potential is the integral of the magnetic fleld, at a fixed 
external magnetic field we have to expect an increasing 
vector potential with increasing penetration depth. The 
increase of the surface vector potential towards low tem- 
peratures has a different physical origin: it is directly 
related to the anomalous Meissner effect and the fleld 
increase shown in Fig. [TOl 

Since the vector potential is proportional to the super- 
fluid velocity, this nonmonotonous temperature depen- 
dence of the vector potential has a direct influence on the 
size of the peak splitting in the local density of statesi^., 
which we show in Fig. [T^l It can be seen that the splitting 



FIG. 11: (Color online) Temperature dependence of the sur- 
face vector potential for orientation a = j. Results are shown 
for K = 10 and Bext = 0.02_Bc2 (triangles), «: = 30 and 
Bext = 0.006-Bc2 (open circles), k = 63 and Bext ~ 0.006-Bc2 
(filled circles). For comparison, the squares show the result 
for K = 10 and Bext = 0.02_Bc2 including impurity scattering 
with a mean free path of I = 5.3^0. Lines are guide to the 
eye. 
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FIG. 12: (Color online) Nonmonotonous splitting of the local 
density of states for orientation a = ^ at different tempera- 
tures. The magnetic field is B = 0.02i3c2. 



is large both for low temperatures and close to Tc- As a 
result also the peak height has a nonmonotonous temper- 
ature dependence. The observation of such an increase of 
the peak splitting towards low temperatures could be an 
experimental signature of the anomalous Meissner cur- 
rents. It should be pointed out, however, that this ef- 
fect becomes less pronounced the larger the k value of 
the material. This is shown in Fig. [TT] for k = 30 and 
K = 63 as the open and solid circles, respectively. For 
these higher values of k the increase of the vector poten- 
tial towards low temperatures is gradually reduced. We 
also want to mention that impurity scattering gradually 
reduces this low temperature increase. The squares in 
Fig. [11] show the behavior for k = 10 and a mean-free 
path of I = 5.3^0- The low temperature increase is re- 
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FIG. 13: Surface vector potential as a function of the external 
magnetic field B/Bc2 for temperature T — O.lTc and k = 10. 
The solid line shows the result for a ~ j and the dashed line 
for a — 0. In the inset the low field range is shown. 
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FIG. 14: (Color online) Local density of states at the surface 
for different impurity mean free paths in the presence of an 
external magnetic field B = 0.006i3c2. The temperature is 
T = O.lTc and the orientation angle q = ^. 



duced, while the increase towards Tc is shifted due to the 
reduction of the bulk Tc- 

The presence of the surface Andreev bound states also 
has a significant influence on the nonlinear Meissner ef- 
fect. We demonstrate this in Fig. [T21 which shows the 
surface vector potential as a function of the external 
magnetic field for orientations a ~ j (solid line) and 
a = (dashed line). For a — 0, where surface Andreev 
bound states are absent, the response is linear over a 
broad range of magnetic fields. In contrast, for a ~ j 
sizeable nonlinear corrections are visible, seen as a steep 
increase at low fields. In the low field range of the order 
of ^ 10^^Bc2, shown in the inset, the response is linear 
in both cases with a significantly larger slope at a = j. 
We suppose that this surface related effect may have an 
important influence on the nonlinear Meissner effect in 
d-wave superconductors and intermodulation distortion 
generated in high- Tc microwave resonators^^iSSii^. 
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FIG. 15: (Color online) Negative imaginary part of the self 
energy E at the surface for different impurity concentrations 
in the presence of an external magnetic field B = 0.006i3c2. 
The temperature is T = O.lTc and the orientation angle a = 




V. INFLUENCE OF IMPURITY SCATTERING 
ON THE PEAK SPLITTING 

Having discussed the two limiting cases of impurity 
scattering without an external magnetic field and the in- 
fluence of a magnetic field without impurity scattering, 
we turn now to a discussion of the combined effect of 
impurity scattering in the presence of an external mag- 
netic field. Naively, one might expect that the impurity 
scattering will wash out the peak splitting. We will show 
below that the situation is more complex, however. In 
the following we will work with a value of k = 63, which 
is more realistic for hole doped high- Tc cuprates. 

In Fig. [14] we demonstrate the infiuence of impurity 
scattering on the peak splitting of the local density of 
states at the surface. Here we have chosen a surface 
angle a = ?, temperature T = O.lTc, and an external 



magnetic field of B = 0.006i?c2- It can be seen that the 
peak height is strongly reduced and the peak width grows 
with decreasing mean free path. However, the size of the 
peak splitting remains almost unaffected. This peculiar 
effect can be understood from the energy dependence of 
the negative imaginary part of the self energy S , which 
is shown in Fig. [151 Here we can see that the scatter- 
ing rate also shows a splitting in energy, which is just 
the mirror image of the splitting in the local density of 
states. It results from the fact that due to the peak split- 
ting the available phase space for scattering processes is 
strongly reduced at low energies. This in turn means that 
the quasiparticle scattering rate at small energies remains 
small even when the mean free path becomes small. This 
leads to a self-stabilization of the peak splitting making 
it robust against impurity scattering. 

In Fig. [injwe show the peak splitting for a fixed mean 
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FIG. 16: (Color online) Local density of states at the surface 
for different external magnetic fields with mean free path I — 
10.6^0. The temperature is T = O.lTc and the orientation 
angle q= f . 
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FIG. 17: (Color online) Imaginary part of the off-diagonal 
self energy with external magnetic field B — 0.006-Bc2 for 
orientation a — j, temperature T = O.lTc and I = 3.5^0. 



free path I = 10.6^o and a series of external magnetic 
fields. When the magnetic field is increased, the peak 
splitting does not evolve as a dip near zero energy, which 
gradually becomes deeper, but instead the splitting opens 
up like a curtain. Note that the frequency dependence 
at low energy is approximately linear over an increasing 
energy scale and the two peaks get a triangular lineshape. 
In the absence of an external magnetic field we have 
pointed out above that the anomalous self energy E^ 



does not vanish except for certain highly symmetric an- 
gles like a = 7r/4. In the presence of an external magnetic 
field even for a — 7r/4, the anomalous self energy 'E^ does 
not vanish anymore. This is shown in Fig. 1171 where the 
imaginary part of S^ is shown as a function of energy 
for different distances from the surface. The reason for 
this is that the special reflection symmetry of the case 
a — 7r/4 is broken now by the direction of the current 
flow. 

VI. CONCLUSIONS 

We have studied the influence of self-consistent Born 
impurity scattering on the surface Andreev bound states 
in a d-wave superconductor. In the absence of an ex- 
ternal magnetic field Born impurity scattering leads to 
a broadening of the zero energy peak in agreement with 
previous worki2iii. We have shown, however, that the 
effect of the impurities is much stronger at the surface 
than in the bulk. Also, the renormalization of the pair 
potential in general does not vanish anymore near the 
surface. 

In the presence of an external magnetic field we have 
demonstrated that for small values of k ~ 10 the zero en- 
ergy peak splitting has a strongly nonmonotonous tem- 
perature dependence. For fixed external magnetic field 
the peak splitting is larger at small and at high temper- 
atures and has a minimum in the intermediate temper- 
ature range. We have also shown that the presence of 
Andreev bound states leads to significant nonlinear cor- 
rections to the Meissner effect. The range of k values 
in question is small compared with hole doped high- Tc 
cuprate materials. However, we suggest that these ef- 
fects may be observable in some electron-doped cuprates, 
which have smaller k values due to a lower Tc and larger 
carrier densities^i. 

We have shown further, that the peak splitting turns 
out to be quite robust against Born impurity scattering. 
This results from a self-stabilizing effect of the peak split- 
ting. 
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